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Abstract 

The reasonableness of the use of perturbative QCD notions in the region close to the scale of 
hadronization, i.e., below < lGeV is under study. First, the interplay between higher orders of 
pQCD expansion and higher twist contributions in the analysis of recent Jefferson Lab (JLab) 
data on the Generalized Bjorken Sum Rule function T{ n (Q 2 ) at 0.1 < Q 2 < 3GeV 2 is studied. 
It is shown that the inclusion of the higher-order pQCD corrections could be absorbed, with good 
numerical accuracy, by change of the normalization of the higher- twist terms. Second, to avoid 
the issue of unphysical singularity (Landau pole at Q = A ~ 400 MeV), we deal with the ghost- 
free Analytic Perturbation Theory (APT) that recently proved to be an intriguing candidate for a 
quantitative description of light quarkonia spectra within the Bethe-Salpeter approach. The values 
of the twist coefficients extracted from the mentioned data by using the APT approach provide 
a better convergence of the higher-twist series than with the common pQCD. As the main result, 
a good quantitative description of the JLab data down to Q ~ 350 MeV is achieved. 

PACS numbers: 11.55.Hx, 11.55.Fv 
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I. INTRODUCTION 



The analysis of Deep Inelastic Scattering (DIS) data by combination of Perturbative 
Quantum Chromodynamics (pQCD) and Operator Product Expansion (OPE) provides us 
with a test site for combining both the perturbative and non-perturbative (NP) QCD con- 
tributions in the low energy domain. In particular, the Generalized (Q 2 -dependent) Bjorken 
Sum Rule (BSR) [lj is a renown target ground for testing different possibilities @, y] • Fortu- 
nately, fresh Jefferson Lab data [4( give information on the spin-dependent BSR amplitude 
T\~ n (Q 2 ) behavior close to the confinement/hadronization scale. Meanwhile, in this region 
the common theoretical pQCD analysis is spoiled by the unphysical singularities in the in- 
frared (IR) region at a scale ~ A ~ 400 MeV . 

To cure this disease (known as the Landau pole trouble) of the pQCD expansion, we use 
the Analytic Perturbation Theory (APT) approach j^] based on the causality principle im- 
plemented as analyticity imperative in the complex Q 2 -plane for the QCD coupling ot s (Q 2 ) 
in the form of the Kallen-Lehmann spectral representation (for a review on APT concepts 
and algorithm see Ref. 0]). 

In principle, the shift of the pQCD frontier (i.e., the boundary above which pQCD is 
applicable) and the rearrangement of the total contribution between perturbative and non- 
perturbative terms is possible by an appropriate modification of perturbative series. Ex- 
amples are provided, say, by IR renormalons 0, ls| and the extractions of higher twists 
(HT) using various approximations of pQCD [9]. In the present paper, we systematically 
explore the possibility of such a shift by extracting the values of HT terms using various 
approximations and modifications of PT to analyse recent JLab data [l(J 11, 12] on the 
BSR. 

We found that a particular form of solution to the Renormalization Group (RG) equation, 



namely, the "denominator" form [13] (see Eq. (j2.4p below) is much more suitable for the 
use in the low-Q 2 region than the most popular ones based on the "multistory" Eq. (9.5) of 
the Particle Data Group [lH compendium (Eq. (12.31) below). The inclusion of higher-order 
(HO) of PT demonstrates a "duality" between HO and HT, in the sense that HT terms 
are absorbed, with good numerical accuracy, into HO terms. As a result, HT coefficients 
decrease. At the same time, we observed that the description of the data is improved only 
up to the two-loop order of PT, which may be a signal of the asymptotic character of PT 
series in the region close to the Landau pole. 

A further shift of the pQCD frontier is achieved by using the analytic APT modifications 
of pQCD in the analysis of BSR below 1 GeV making possible investigation of the interplay 
of PT and HT contributions. As a result, we find that while the next-to-leading twist 
(114) term is larger in APT than in the usual pQCD, the HT (fJ,e,s) coefficients in APT are 
smaller (making the impression of HT series being convergent) and allowing for a reasonable 
description of the data down to Q ~ 350 MeV. 



II. THE BJORKEN SUM RULE IN CONVENTIONAL PT 

The Bjorken integral is defined 

rr n (Q 2 ) = C dxMx, Q 2 ) - fi(x, Q 2 )) , (2.1) 
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via the spin- dependent proton and neutron structure functions , with x = Q 2 /2Mu 1 
the energy transfer v and the nucleon mass M. At large Q 2 , the BSR comes to its renowned 
form T\~ n = qaI^i where qa = 1.267 ± 0.004 is the nucleon axial charge defined from the 
neutron /3-decay. At finite Q ^> A , the BSR is given by the OPE series in 1/Q l ~ 2 with even 
i = 2,4... being the number of a twist and the pQCD series in a™ . The expression for the 
perturbative part of T\~ n (Q 2 ) including the HT contribution is (see e.g. Ref. 
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with numerical values given at rif = 3 and weak dependence of [i 2 i on logQ 2 neglected. The 
first non-leading twist term 15j can be expressed [12] 



,,P- n ~ 4Af 2 fP-n 



in terms of the colour polarizability f 2 ■ 



Within the pQCD, the a s coupling is usually taken in the form (Eq. (7) in [10], Eq. (9.5) 



in PDG [l4|) expanded over In L/L with (L = ln(Q 2 /A 2 ), b k = p k /ft 
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Here, the second term in the first line is the 2-loop contribution and the framed term 
usually is referred to as "the 3-loop one" , while all contents of the second line is treated on 
the equal footing with the 4-loop term ~ b 3 . 

However, it is evident that pieces of genuine 2-loop contribution proportional to b\ are 
entangled with the higher-loop ones. This defect is absent in the more compact "Denominator 
representation' 
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which, being generic for the PDG expression, is closer to the iterative RG solution and, 
hence, more precise. Below, we shall refer to it as to "Denom" . 

A detailed higher-twist analysis based on the total set of low energy SLAC and JLab data 
was performed in Ref. [ll[. The result of the combined fit done in the Q 2 -range 0.66-10.0 
GeV 2 is f 2 (Q 2 = lGeV 2 ) = -0.101 ± 0.027 and /WM 4 = 0.084 ± 0.011 (elastic contribution 
included). The fitting procedure of Refs. [HI, I3| taking the pQCD leading- twist term 
calculated at NLO a s and using the 2-loop "Denom" coupling a^' D was repeated. We 
succeeded in obtaining the central values of Refs. 11 , l3] in the two-parametric fit with the 
output 

f 2 = -0.096 ± 0.012, /x 6 /M 4 = 0.087 ±0.004, x 2 = 0.48, 

where the errors are statistical only. These results are compatible with HT extraction 
performed in Ref. 17 . 
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It is of special interest to study the BSR data with the elastic contribution (necessarily 
present in the OPE framework [l8() excluded, since the low-Q 2 behavior of such an "in- 
elastic" BSR integral (coinciding with the usual BSR for Q 2 — > oo) is constrained by the 



Gerasimov-Drell-Hearn (GDH) sum rule [19j, and one may investigate its continuation to low 
energies j3|. To this goal, doing the same NLO combined fit (elastic contribution excluded), 
one gets 



inel 



-0.080 ± 0.016, /V M = 0.022 ± 0.005, x = 0.91 



The difference is noticeable starting from which is natural due to a decrease of an elastic 
contribution with growing Q 2 . 

To explore the fit results sensitivity to the PT order and to the form of a s below 1 GeV, 
we have also performed fits at the 1-, 2- and 3- loop levels. The minimal borders of fitting 
domains in Q 2 were settled from the ad hoc restriction \ 2 ^ 1. 

From Fig. [H one sees that the results obtained with 2- and 4-loop expressions for the 
"Denom" coupling are better consistent with the BSR data at Q < 1 GeV than those based 
on the PDG expression (12.31) . though 3- loop results do not differ significantly. 



0,22 
0,20 ■ 
0,18 
0,16 
0,14 
0,12 
0,10 
0,08 
0,06 
0,04 
0,02 


-0,02 



p-n 



N LO "Denom" 

N 3 LO PDG 

N 2 LO "Denom" 

N 2 LO PDG 

NLO "Denom" 

NLO PDG 

LO Bjorken 



JLab Hall B (CLAS EG1b) [2006] 
(stat. errors only) 

JLab Halls A,B E94010/EG1a [2002] 
JLab Hall B (CLAS EG1a) [2003] 
SLAC E143 [2000] 




0,9 1,0 

Q, GeV 



FIG. 1: Best 3-parametric fits of JLab and SLAC data on Bjorken SR calculated within a s in PDG form 
(|2.3p and in the "Denom" one (|2.4p at various loop orders. 

Indeed, by using Eq. (12.41) one may extend the applicability domain of Eq. (12. 2p down to 
Q 2 ~ 0.27 GeV 2 . At the same time, Eq. (TOD works well only down to ~ 0.47 GeV 2 due to 
extra In™ L/L n+l singularities. 

The fitting of BSR data in 2,3,4-loops over the fixed range 0.6 GeV < Q < 2.0 GeV yields 
a "swap" between the higher orders of PT and HT terms. In Fig. [21 we show one-parametric 
fits with 2,3-loop a s pQCD to the BSR. One can see there that the higher-loop contributions 
are effectively "absorbed" into the value of /X4 which magnitude decreases with increasing 
of the loop order while all the fitting curves are very close to each other. This observation 
reveals a kind of "duality" between perturbative a s -series and nonperturbative l/Q 2 -series. 

This also means the appearance of a new aspect of quark hadron duality, the latter being 
the necessary ingredient of all the QCD applications in the low energy domain. Usually, it 
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is assumed 2(| that the perturbative effects are less important there than the power ones 
due to a nontrivial structure in QCD vacuum. 

In our case, the PT corrections essentially enter into the game, so that the pQCD HO 
terms are relevant in the domain where the concepts of traditional hadronic physics are 
usually applied. 
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FIG. 2: One-parametric fits of the JLab and SLAC data on Bjorkcn SR calculated in the "Denom" form 
in different loop orders. 



The interplay between partonic and hadronic degrees of freedom in the description of 
GDH SR and BSR may be also observed in the surprising similarity between the results of 
"resonance" (Hf and "parton" [3| approaches. 

At the same time, from Fig. [1] it follows that the higher (3- and 4-loop) PT orders yield a 
worse description of the BSR data, probably implying the asymptotic character of the series 
in powers of a s . 

One may ask to what extent the troubles mentioned above are due to the unphysical 
singularities at Q ~ A in PT series for T^~p T . Their influence is essential just at Q < 1 GeV 
where the HT terms play an important role. 

The APT is free of such problems, thus providing a tool to investigate the behavior of 
HT terms extracted directly from the low energy data. This provides a motivation for the 
analysis performed in the next section. 
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III. THE BJORKEN SUM RULE IN APT 



According to the approach developed by Igor Solovtsov and co-authors 22|, the APT 
modification of BSR with HT power corrections looks like 
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It should be noted that the APT Euclidean functions in the 1-loop case are simple enough 
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the higher .4^ being related to the lower ones recursively by differentiating. Analogous 2- 
and 3- loop level expressions involve a little known special Lambert function and are more 
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intricate 

Meanwhile, even for the 3-loop APT case, there exists a possibility to employ the effective 
log approach proposed by Igor Solovtsov and one of the authors [25J and extended recently 
(see Eqs. (12), (14) in Ref. " 
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to higher APT Euclidean and Minkowskian functions. In 
the present context, one may use simple model one- loop expressions (13.21) with some effective 
2-loop log L* accumulating the 2-loop log-of-log 
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L* =L + B(n f )hx Vl 2 + 2vr 2 ; B = (3 1 /f3 2 . 



Happily enough, the second term does not undergo a significant variation in the intermediate 
few GeV region. Indeed, as B{rif = 3) = 0.79 and -8(4) = 0.74, in the region Q < 5 GeV 
involved in the BSR data analysis, a simple approximation 
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happens to be accurate enough. That is, instead of the cumbersome 3-loop expressions for 
the APT functions, in Eq. (13. ip one can use the 1-loop expressions (13.21) with A mo d 
value given by the second relation (13 .4p . 

If we take A (3) = 380 MeV j27[, then A^ ~ 190 MeV. The corresponding maximal 
errors of the model (13.31) for first and second functions are j26| 8A mod /A™ od — 4% and 
5A n f od /A^ od ~ 8% at Q ~ A® , which seems to be sufficiently accurate. Indeed, as far as 
Ai{Q = 400 MeV) = 0.532 and A 2 (400MeV) = 0.118, the total error in r^ n PT is mainly 
determined by the first term, being of the order ST p ~ n /T p ~ n ~ 5A 1 f od /'K ~ 1%, i.e., less 
than the data uncertainty. 
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Turn now to the 3-loop APT part of the Bjorken integral r^~A PT (Q 2 ). Its value is quite 
stable with respect to small variations of A (in contrast with huge instability of T^~p T ): it 
changes now by about 2 - 3% within the interval A (3) = 300 - 400 MeV 1 . 

We also performed the comparison with Simonov's "glueball-freezing model" (SGF- 
model) - see Fig. [3], with similar to PDG 1/L-type loop expansion for freezed coupling 

a B (Q 2 ) = af\L) , L = ln(^J^) (SGF) 

where 2-loop a s is taken in the form of two terms from the first line in Eq. (12.31) with 
logarithm modified by a "glueball mass" Mo ~ 1 GeV and the usual PT expansion in powers 
of as in T p ~ n is adopted. 

Extending the analysis of Ref. [22| to lower Q values, we estimated the relative size of 
APT contributions to the BSR. It turned out that the third term ~ ^3 contributes no more 
than 5% to the sum, thus supporting the practical convergence of the APT series. 

Note that the APT functions Ak contain the (Q 2 )~ k power contributions which effectively 
change the fitted values of //-terms. In particular, subtracting extra (Q 2 ) -1 term induced 
by the APT series 



with x = 0.43, we get 



^if" ^-0-048, A«~0.2GeV, (3.5) 



that nicely correlates with the result in Ref. llj: /14/M 2 ~ —0.045. This demonstrates the 
concert of the APT analysis with the usual PT one for the BSR data at Q 2 > 1 GeV 2 . 

In Fig. [31 we show best fits of the combined data set for the function T\~ n {Q 2 ) in the PT 
and the APT approaches. The corresponding numerical results are given in Table [H Our 
fit gives the HT values indicating a better convergence of the OPE series due to decreasing 
magnitudes and alternating signs of consecutive terms, in contrast to the usual PT fit results. 



It is worth noting that the best APT fit allows one to describe well all the BSR data 
at scales down to Q ~ 350 MeV with only the first three terms of the OPE series, unlike 
the usual PT case, where such fits happened to be impossible (due to the ghost issue) even 
for an increasing number of HT terms. This means, that the lower bound of the pQCD 
applicability (supported by power HT terms) now may be shifted down to Q ~ Aq CD ~ 350 
MeV. 

However, it seems to be difficult to get a description in the region Q < Aqqp. This is 
not surprising, because the expansion in positive powers of Q 2 and its matching [3j with the 
HT expansion are relevant here. In this respect, Aqcd scale appears as a natural border 
between "higher twist" and "crural" nonperturbative physics. 



In particular, this means that the low-Q BSR data cannot be used for a determination of A in the APT 
approach. 
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FIG. 3: Best 1,2,3-paramctric fits of the JLab and SLAC data on Bjorkcn SR calculated with NLO "Denom" 
(solid lines) and PDG (dashed lines) couplings and N 2 LO APT (dash-dotted lines) at fixed Aqcd value 
corresponding to the world average. We also show the pQCD part of the BSR at different values of A^ 3 -* 
300, 350, 400 MeV calculated within APT (short-dashed lines) and SGF-model 
the glueball mass M = 1.2, 1.0, 0.8 GeV (with A = 380 MeV) (dotted lines). 



at different values of 



TABLE I: Combined fit results for the HT terms in APT and conventional PT in PDG and "Denomi- 
nator" forms. 



Method 


QLm GeV 2 


/x 4 /M 2 


/i 6 /M 4 


/, 8 /M 6 


NLO PDG 


0.50 


-0.043(2) 








A = 380 MeV 


0.30 


-0.074(4) 


0.025(2) 







0.27 


-0.049(5) 


-0.007(5) 


0.009(1) 


NLO "Denom" 


0.47 


-0.046(2) 








A = 340 MeV 


0.17 


-0.066(2) 


0.013(4) 







0.17 


-0.061(4) 


0.009(3) 


0.0005(3) 


N 2 LO APT 


0.47 


-0.054(1) 








A = 380 MeV 


0.17 


-0.065(2) 


0.0081(5) 







0.10 


-0.069(2) 


0.0114(9) 


-0.0006(1) 



IV. CONCLUSION AND OUTLOOK 

The separation of perturbative and NP physics may be different if some modification of 
perturbation theory is adopted. To test such a separation, we performed a systematic com- 
parison of HT terms extracted from the very accurate JLab data on BSR in the framework 
of both the common PT and yhe APT in QCD and came to the following results. 
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• The evidence that the "Denominator" form (12. 4p of the QCD coupling a s is more 
suitable in the low Q region is given (see Fig. 1). 

• A kind of duality between HO of PT and HT is observed so that HO terms absorb part 
of HT contributions moving the pQCD frontier between the PT and HT contribution 
to lower Q values (see Fig. 2). 

• This situation is more pronounced in the APT where convergence of both the HO and 
HT series is much better. While the twist-4 term happened to be larger in magnitude 
in the APT than in the PT, the subsequent terms are essentially smaller and quickly 
decreasing (as the APT absorbs some part of NP dynamics described by HT) . This is 
the second reason of the shift of pQCD frontier to lower Q values. 

As the main result, a satisfactory description of the data down to Q ~ ^qcd — 
350 MeV is achieved by taking the analytic HO and HT contributions into account 
simultaneously (see Fig. 3). 

In a sense, this could be natural if the main reason of such a success was the disappearance 
of unphysical singularities. We have in mind that the singularity-free APT and SGF QCD 
couplings are ver y c lose in the domain Q > 400 MeV. Moreover, various lattice data [2£| 
(see also reviews 3(| and references therein) yield similar a s curves there. 



It will be very interesting to explore the interplay between PT and NP physics against 
other low energy experimental data. 
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